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Abstract 

In the CP-invariant super symmetric theories, the steep S-wave (slow P-wave) rise 
of the cross section for any non-diagonal neutralino pair production in e + e~ anni- 
hilation, e + e~ — » XiX^j {i j)i uear threshold is accompanied by the slow P-wave 
(steep S-wave) decrease of the fermion invariant mass distribution of the 3-body 
neutralino decay, Xi ~~ Xjff (f = ' or Q)i near the end point. These selection 
rules, unique to the neutralino system due to its Majorana nature, guarantee that 
the observation of simultaneous sharp S-wave excitations of the production cross 
section near threshold and the lepton or quark invariant mass distribution near the 
end point is a qualitative, unambiguous evidence for CP violation in the neutralino 
system. 
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1 Introduction 



Most supersymmetric extensions of the Standard Model (SM) based on some soft super- 
symmetry (SUSY) breaking mechanism contain several CP phases, whose large values 
tend to render lepton and quark electric dipole moments (EDM) too large to satisfy strin- 
gent experimental constraints pp. Such CP crises are generic in supersymmetric theories, 
but may be resolved by pushing the masses of some sparticles, especially the first and 
second generation sfermions, above a few TeV, by arranging for internal cancellations, or 
by simply setting phases to be extremely small j2]. On the other hand, new sources of 
CP violation beyond the SM are required to explain the non-zero baryon asymmetry in 
the universe in the standard Big Bang framework j3J. Therefore, it is crucial to look for 
new signatures for CP violation in such SUSY scenarios with some large phases, as long 
as they are consistent with the stringent EDM and other low-energy constraints. In this 
light, detailed analyses of the neutralino sector at future e + e~ linear collider experiments 
jl] can prove particularly fruitful HI El M > because in most supersymmetric theories 
neutralinos belong to the class of the lighter supersymmetric particles ^0] arid the neu- 
tralino system contains two non-trivial CP violating phases. 

There are many different ways for probing CP violation in the neutralino system. The 
imaginary parts of the complex parameters in the neutralino mass matrix could most 
directly and unambiguously be determined by measuring suitable CP violating observ- 
ables by exploiting initial beam polarization and angular correlations between neutralino 
production and decay at future high-energy colliders 00 El HI- But, their experimental 
measurements will be quite difficult. The presence of the CP violating phases can also 
be identified through by their impact on CP-even quantities such as neutralino masses, 
branching ratios and so on. However, since these quantities are already non-zero in the 
CP conserving case, the detection of the presence of non-trivial CP phases will require a 
careful quantitative analysis of a number of physical observables, especially for small CP- 
odd phases giving rise to very small deviations from the CP-conserving values [I]. On the 
other hand, the rise of excitation curves near threshold for non-diagonal neutralino pair 
production in e + e~ collision is altered qualitatively in CP-noninvariant theories [SHE], by 
allowing the steep S-wave increase of all pairs simultaneously. Thus, as demonstrated in 
Ref. precise measurements of the threshold behavior of the non-diagonal neutralino 
pair production processes may give clear indications of non-zero CP violating phases in 
the neutralino sector, if at least three different neutralino states are accessible kinemati- 
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cally. 

In the present note we provide a new powerful method for probing CP violation in 
the neutralino system, which is based on a combined analysis of the threshold excita- 
tions of neutralino pair production in e + e~ annihilation and the fermion invariant mass 
distribution near the end point of the 3-body neutralino fermionic decays: 

e+e-^xlx] and %^%ff (J = l,q). 

[The 3-body decay process includes clean and e + e~ decay channels with little 

background, which allow a clear reconstruction of the kinematical configuration with 
good precision.] This method relies on selection rules, unique to the neutralino system 
due to its Majorana nature in CP-invariant theories, and it can work effectively if the 
branching ratios of the 3-body neutralino fermionic decays are not suppressed. [Once 
two-body decays of the neutralino xl Higgs bosons or sfermions are open, the 

new method is ineffective.] 

Before demonstrating the new method for probing CP violation in the neutralino 
system in detail, we describe briefly the mixing for the neutral gauginos and higgsinos 
in CP-noninvariant theories with non-vanishing phases in Sec. 121 In Sec. El we introduce 
the selection rules for the production of neutralino pairs and the neutralino to neutralino 
transition via a (virtual) vector boson or sfermion exchange. Then, we prove that in 
any CP-invariant SUSY theory, if the production cross section for any non-diagonal 
neutralino pair in e + e~ annihilation increases steeply in S-waves (slowly in P-waves) 
near threshold, the lepton or quark invariant mass distribution of the decay Xi ~ * X°j ff 
(/ = / or q) decreases slowly in P-waves (steeply in S-waves) near the end point. Thus, the 
observation of simultaneous sharp S-wave excitations of both the production of any non- 
diagonal neutralino pair XiXj near threshold and the fermion invariant mass distribution 
of the decay x° — > x° // near the end point will be a qualitative, unambiguous evidence 
for CP violation in the neutralino system. A quantitative demonstration of the method 
based on a specific set of the relevant supersymmetry parameters is given in the last part 
of Sec. 01 Finally, conclusions are drawn in Sec. |U 

2 Neutralino Mixing 

In the minimal supersymmetric extension of the Standard Model (MSSM), the mass 
matrix of the spin-1/2 partners of the neutral gauge bosons, B and W 3 , and of the 
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neutral Higgs bosons, H® and H 2 , takes the form 
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in the {B, W 3 , H®, if^ } basis. Here Mi and M 2 are the fundamental supersymmetry 
breaking U(l) and SU(2) gaugino mass parameters, and /i is the higgsino mass parameter. 
As a result of electroweak symmetry breaking by the vacuum expectation values of the two 
neutral Higgs fields V\ and v 2 {sp = sin/5, cp = cos/3 where tan/3 = v 2 /vi), non-diagonal 
terms proportional to the Z-boson mass m z appear and the gauginos and higgsinos mix 
to form the four neutralino mass eigenstates Xi {i — 1-4). In general the mass parameters 
Mi, M 2 and fi in the neutralino mass matrix can be complex. By re-parameterization 
of the fields, M 2 can be taken real and positive, while the U(l) mass parameter Mi is 
assigned the phase $i and the higgsino mass parameter /i the phase 

The neutralino mass eigenvalues m ; = m^o (i = 1,2,3,4) can be chosen positive 
by a suitable definition of the mixing matrix N, rotating the gauge eigenstate basis 



{B, W 3 , H®, H 2 } to the mass eigenstate basis of the Majorana fields x 



1-4). In 



general the matrix iV involves 6 angles and 10 phases, and can be written as [HI E] 

N = diag {e ia \ e ia \ e ias , e ia * } R 34 R 24 Ri4 R 23 Ris R12 , (2) 

where Rjk are rotations in the complex [jk] plane characterized by a mixing angle 9jk and 
a (Dirac) phase (3jk- One of (Majorana) phases is nonphysical and, for example, a\ may 
be chosen to vanish. None of the remaining 9 phases can be removed by rotating the fields 
since neutralinos are Majorana fermions. The neutralino sector is CP conserving if [i and 
Mi are real, which is equivalent to j3ij = (mod tt) and ai = (mod vr/2). Majorana 
phases a« = ±7r/2 do not signal CP violation but merely indicate different intrinsic CP 
parities of the neutralino states in CP-invariant theories 



3 Neutralino Pair Production and 3— Body Decays 

Both the production processes, e + e~ X^x'j [hj = 1~4), and the 3-body neutralino 
decays, Xi ~^ Xj J J "> are generated by the five mechanisms: s-channel Z exchange, and t— 
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and w-channel fi,R exchanges with / = e for the production processes. After appropriate 
Fierz transformations of the sfermion exchange amplitudes and with the fermion masses 
neglected, the transition matrix element of the production process e + e~ — > XiXj an d that 
of the 3-body fermionic neutralino decays x° — ► X° // can be written as 

T(e + e-^tx°) = E Q*p[He + h^au{e-)][v($WPpv(j$], (3) 

a,P=L,R 

D(Xi^X°jff) = E Q'ap [uifh^vif)] [uix^Ppuixl)] , (4) 

a,P=L,R 

that is to say, as a sum of the products of a x° vector or axial vector current and a 
fermion vector or axial vector current, respectively. We refer to Ref. [Hj and Ref. [Hj for 
the expressions of the generalized bilinear charges Q a/ 3 and Q' a/3 , just mentioning that the 
bilinear charges become independent of the kinematical variables when two neutralinos 
are at rest. Therefore, in this static limit, both the production and the decays can be 
considered to proceed via a static vector boson exchange. 

Some general properties of the bilinear charges Q a p and Q' a/3 in Eqs. dHJ) and (jlj can be 
derived in CP-invariant theories by applying CP invariance and the Majorana condition 
for neutralinos to the transition matrix elements. In CP-invariant theories, the production 
of a neutralino pair through a vector or axial vector current with positive intrinsic CP 
parity satisfies the CP relation US] 

l = rM-l) L , (5) 

in the non-relativistic limit of two neutralinos, where rj l = ± i is the intrinsic CP parity 
of Xi an d L is the orbital angular momentum of the neutralino pair. The selection rule 
((HJ) reflects the fact that if two neutralinos Xi an d xl have the same or opposite CP 
parity, the current for the neutralino pair production must be pure axial-vector or pure 
vector form, respectively, cf . ^3] • Because the axial-vector current and the vector current 
involve the combination of u and v spinors for the two Majorana particles, the axial vector 
corresponds to the P-wave (L = 1) and the vector to the S-wave (L — 0). 

On the other hand, the neutralino decay, Xi ~ > x'j + V, where V stands for the final 
fermion current in Eq. satisfies the CP relation 

rf = r] j (-l) L or equivalently 1 = -rfrf (-1) L , (6) 

in the non-relativistic limit of two neutralinos, where L is the orbital angular momentum 
of the final state of an d V. We emphasize first that the neutralino to neutralino 
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transition current is pure axial-vector or pure vector form for the two neutralinos of the 
same or opposite CP parity, respectively, as in the production case. However, because 
two M-spinors are associated with the currents in the neutralino to neutralino transition, 
the axial-vector corresponds to S-wave excitation while the vector corresponds to P-wave 
excitation, giving rise to the relative minus sign between (JSJ) and (jHJ. 

One immediate consequence of the selection rules (0) and © is that, in CP-invariant 
theories, if the production of a pair of neutralinos with the same (opposite) CP parity 
through a vector or axial vector current is excited slowly in P-waves (steeply in S-waves) 
, then the neutralino to neutralino transition via such a vector or axial vector current 
is excited sharply in S-waves (slowly in P-waves). More explicitly, the power of the 
selection rules (JHJ) and (JOJ) can clearly be seen by inspecting the expressions for the S- 
wave excitations of the total cross section <r{ij} (i 7^ j) near threshold and of the fermion 
invariant mass distribution of the 3-body neutralino decay Xi ~ y x'j ff (with the fermion 
masses neglected) near the end point: 

^ Wm t m r ^ 2 + ^ j + 

(mi + mjf I J 

dT{ij} _ 2a 2 {nij\ 3 / 2 . ... , ? , 1 , ,. r „ ,2 , . w -,, ,2 1 „ v *» 



dzff n \rrii 



m { - mj) ft' j I ^eG' R \ 2 + | 3teG'J 2 } + 0(ft K 



where ft = y 1 — (mj + mj) 2 / s and ft' = Jl — z 2 j with the dimensionless variable z// = 
m /// m / 1 / x > the ratio of the fermion invariant mass mjf to its maximal value m™j x = 
mi — my Here, the coupling dependent parts, each of which is connected with the chirality 
of the neutralino current, are given by 



G { £ = {I \ 2 Qf f w) D(>\N i3 N* 3 - N i4 N* A ) + -^if iV^* , (9) 
^ C W S W s w c w 



with N' i2 = (ig — Qf) swNn — j| cwN i2 , and the kinematic functions are given by 



D = (mi + m 3 -) 2 / ((mj + m^) 2 - m|) , 
Fl,r = {mi + mj) 2 / {m\ h R + mmi^j) , 
D' = (mi - mj) 2 / ((mi - mj) 2 - m 2 z ) , 

f l, r = (mi - mj) 2 / {mj - m^) . (10) 
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In CP-invariant theories, all the (complex) rotation matrices Rjk in Eq. (J2J) become real 
and orthogonal. Therefore, if the neutralinos Xi an d X°j have the same CP parity, then the 
Major ana phase difference, or 7r, and so N^N^ is real. On the contrary, if the 

neutralino pair have the opposite CP parity, the phase difference ±7r/2 and so 

NikN* L is purely imaginary. Consequently, in CP-invariant theories the cross section of a 
non-diagonal neutralino pair rises steeply in S-waves only when the produced neutralinos 
have the opposite parity, as dictated by the first CP relation (JHJ) and as clearly indicated 
by Eq. (|7jl. One important implication of the selection rule is that, even if the {ij} and 
{ik} pairs are excited steeply in S-waves, the pair {jk} must be excited slowly in P-waves 
characterized by the slow rise ~ /3 3 of the cross section [13 El- In contrast to the production 
case, the characteristic sharp S-wave decrease of the fermion invariant mass distribution 
near the end point is possible only if the neutralinos have the same CP parity, as dictated 
by the second CP relation (JHJ) and as clearly indicated by Eq. (jHj). 

However, in the CP-noninvariant theories the orbital angular momentum is no longer 
restricted by the selection rules (0) and (JHJ). The production of all non-diagonal pairs 
can simultaneously be excited steeply in S-waves near threshold, and the corresponding 
neutralino to neutralino transition can be excited steeply in S-waves even if the production 
cross section of the same non-diagonal neutralino pair is excited steeply in S-waves. 
Consequently, CP violation in the neutralino system can clearly be signalled by (i) the 
sharp S-wave excitations of the production of three non-diagonal {ij}, {ik} and {jk} pairs 
near threshold [H] or by (ii) the simultaneous S-wave excitations of the production of any 
non-diagonal {ij} neutralino pair in e + e~ annihilation, e + e~ — > near threshold and 

of the fermion invariant mass distribution of the neutralino 3-body decays, x°i ~~ * //> 
near the end point. 

It is noteworthy that only the light neutralinos Xi 2 among the four neutralino states, 
which are expected to be lighter than sfermions and gluino in many scenarios, may be 
kinematically accessible in the initial phase of e + e~ linear colliders. In this situation, 
the method based on the threshold behaviors of the production of three different non- 
diagonal neutralino pairs for probing CP violation is not available. On the contrary, the 
combined analysis of the threshold excitation of the production process, e + e~ — > X1X21 
and the fermion invariant mass distribution of the decay, X2 ~~ > Xiffi near the end point 
can still serve as one of the most powerful probes of CP violation in the neutralino system 
even in the initial phase of e + e~ linear colliders. 
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Figure 1: (a) The threshold behavior of the neutralino production cross-sections c{12] 
near the threshold and (b) the lepton invariant mass distribution of the decay x° ~~ > Xi 
near the end point, illustrated for the parameter set: tan/5 = 10, \Mi\ = 100 GeV, 
M 2 = 150 GeV, \fi\ = 400 GeV and $^ = as well as the slepton masses, mj L = 250 
GeV and m r = 200 GeV. 



In order to illustrate the method for probing CP violation numerically, we take a 
parameter set for the fundamental SUSY parameters 1 : 

tan/? =10; |Mi| = 100 GeV, M 2 = 150GeV, |/i|= 400GeV; $ M = (11) 

and we choose two different values, {0, 7r} for the phase $i, in the CP-invariant case and 

one value, 7r/2, in the CP non-invariant case. [The parameter point with such a large 

phase $i = 7r/2 might already have been excluded by the stringent EDM constraints. 

Nevertheless, this point is taken just for illustrative purpose in the present work; the 

indirect EDM limits depend also on many parameters of the theory outside the neutralino 

sector.] We take the slepton masses, mj L = 250 GeV and mj R = 200 GeV and consider the 

3-body leptonic decay x° "~ > Xi especially with I = e, fi, for the illustration. We note 

that the neutralinos Xi an d X2 have the same (opposite) CP parity for $1 = ($1 = 71). 

As expected from the selection rules and © in the CP-invariant case, Figure ^ clearly 

shows that if the production of the neutralino pair X1X2 * n e+e ~ annihilation increases 

slowly in P-waves (steeply in S-waves) near threshold, then the lepton invariant mass 

distribution of the decay X2 ~^ Xi decreases steeply in S-waves (slowly in P-waves) 

1 Analyses of electric dipole moments strongly suggest that CP violation in the higgsino sector will be 
very small in the MSSM if this sector is non-invariant at all 
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near the end point for the neutralino pair of the same (opposite) CP parity with $1 = 
($1 = 7i"). On the contrary, in the CP-noninvariant case (<£>i = 7r/2) the production and 
decay are excited steeply both in S-waves. 



4 Conclusions 

We have shown that only in CP-noninvariant theories the production of any non-diagonal 
neutralino pair XiX'j 7^ j) i n e + e~ annihilation near threshold and the fermion invariant 
mass distribution of the 3-body neutralino fermionic decay Xi ~ y X°j ff near the end point 
can simultaneously be excited steeply in S-waves. 

In light of the possibility that only the two light neutralinos Xi an d X2 among the four 
neutralinos can be accessed kinematically in the initial phase of e + e~ linear colliders, the 
combined analysis of the production of the neutralino pair X1X2 i n e + e~ annihilation near 
threshold and the neutralino decay x° ~~ > Xi ff near the end point of its fermion invariant 
mass could provide a first qualitative indication of the CP violation in the neutralino 
system. 

Acknowledgments 

The author would like to thank M. Drees, J. Kalinowski and P.M. Zerwas for helpful 
comments. The work was supported in part by the Korea Research Foundation Grant 
(KRF-2002-070-C00022) and in part by KOSEF through CHEP at Kyungpook National 
University. 



References 

[1] S. Pokorski, J. Rosiek and C.A. Savoy, Nucl. Phys. B570 (2000) 81 piep-ph/9906206| ; 
V. Barger, T. Falk, T. Han, J. Jiang, T. Li and T.Plehn, Phys. Rev. D 64 (2001) 



056007 |hep-ph/ 0101106|. 



[2] Y. Kizukuri and N. Oshimo, Phys. Lett. B249 (1990) 449; T. Ibrahim and P. Nath, 
Phys. Rev. D 58 (1998) 111301 [Erratum-ibid. D 60 (1998) 099902] |hep-ph/9807"50l| 
and hep-ph/0107325l T. Ibrahim and P. Nath, Phys. Rev. D 61, 093004 (2000) 



hep-ph/9910553|; M. Brhlik, G.J. Good and G.L. Kane, Phys. Rev. D 59 (1999) 



9 



115004; R. Arnowitt, B. Dutta and Y. Santoso, Phys. Rev. D 64 (2001) 113010 
hep-ph/0106089|; T. Ibrahim and R Nath, |arXiv:hep-ph/0210"25"l"l 

[3] A.G. Cohen, D.B. Kaplan and A.E. Nelson, Ann. Rev. Nucl. Part. Sci. 43 (1993) 27; 
B.M. Gavela, P. Hernandez, J. Orloff, O. Pene and C. Quimbay, Nucl. Phys. B430 
(1994) 382; V.A. Rubakov and M.E. Shaposhnikov, Usp. Fiz. Nauk 166 (1996) 493. 

[4] TESLA Technical Design Report, Part: III Physics at an e + e~ Linear Col- 
lider, eds. R.-D. Heuer, D. Miller, F. Richard and P. Zerwas, DESY 2001-011, 
arXiv:hep-ph/0106315; T. Abe et al. [American Linear Collider Working Group Col- 
laboration], "Linear collider physics resource book for Snowmass 2001. 2: Higgs and 
supersymmetry studies," in Proc. of the APS/DPF/DPB Summer Study on the Fu- 
ture of Particle Physics (Snowmass 2001) , ed. N. Graf, |arXiv:hep- ex/0106056 K. 
Abe et al, JLC Roadmap Report, presented at the ACFA LC Symposium, Tsukuba, 



Japan 2003, http://lcdev.kek.jp/RMdraft/ 



[5] S.T. Petcov, Phys. Lett. B178 (1986) 57. 

[6] SY. Choi, J. Kalinowski, G. Moortgat-Pick and P.M. Zerwas, Eur. Phys. J. C 22 
563 (2001) |hep-ph/0108117| , ibid. C 23 769 (2002) |hep-ph/0202039"| ; J. Kalinowski, 
|arXiv:hep-ph/0306272l 

[7] G. Moortgat-Pick and H. Fraas, Phys. Rev. D 59 (1999) 015016 |hep-ph/9708"48l] ; 

G. Moortgat-Pick, H. Fraas, A. Bartl and W. Majerotto, Eur. Phys. J. C 9 (1999) 
521 [Erratum-ibid. C 9 (1999) 549] |hep-ph/9903"220| ; G. Moortgat-Pick, A. Bartl, 

H. Fraas and W. Majerotto, Eur. Phys. J. C 18 (2000) 379 |hep-ph/0007222| . 

[8] SY. Choi, H.S. Song and W.Y. Song, Phys. Rev. D 61 (2000) 075004 
|hep-ph/9907474| . 

[9] G. Moortgat-Pick and H. Fraas, Eur. Phys. J. C 25 (2002) 189 |hep-ph/0204"333| . 

[10] See, for instance, B.C. Allanach et al, Eur. Phys. J. C 25 (2002) 113 
|hep-ph/0202233| . 

[11] L. Chau and W. Keung, Phys. Rev. Lett. 53 (1984) 1802; M. Gronau and J.Schechter, 
Phys. Rev. Lett. 54 (1985) 385 [Erratum-ibid. 54 (1985) 1209]; H. Fritzsch and J. 
Plankl, Phys. Rev. D 35 (1987) 1732; F. del Aguila and J.A. Aguilar-Saavedra, Phys. 



10 



Lett. B386 (1996) 241 |hep-ph /9605418|; J.A. Aguilar-Saavedra and G.C. Branco, 
Phys. Rev. D 62 (2000) 096009 lhep-ph/0007025] . 

[12] J. Ellis, J.M. Frere, J.S. Hagelin, G.L. Kane and S.T. Petcov, Phys. Lett. B 132 
(1983) 436. 

[13] B. Kayser, Phys. Rev. D 26 (1982) 1662; B. Kayser and A.S. Goldhaber, Phys. 
Rev. D 28 (1983) 2341; B. Kayser, Phys. Rev. D 30 (1984) 1023; B. Kayser, F. 
Gibrat-Debu and F. Perrier, World Sci. Lect. Notes Phys. 25 (1989) 1. 



11 



